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^^ I The predictions of the IBM two-body random ensemble are compared 

to empirical results on nuclei from Z=8 to 100. Heretofore unrecognized but 
robust empirical trends are identified and related both to the distribution 
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D ■ of valence nucleon numbers and to the need for and applicability of specific, 

^^ , non-random interactions. Applications to expected trends in exotic nuclei are 

^ i discussed. 
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Recently there has been renewed interest in the spectroscopic properties of Hamiltonians 



t^^ ! with random interactions ||T|-0]. Scrutiny has focused on nuclear systems but applications 

O I exist in other many body quantal systems such as quantum dots or metallic clusters. In the 

(^ ' nuclear case, analyses have been carried out in the frameworks of the Shell Model |jll-§,|^ 

O ■ and the Interacting Boson Model (IBM) 0-§[. 

These studies are extremely important because they explore the origins of fundamental 

features of nuclear structure. The upshot of these studies has been surprising recognition 

that some of the most hallowed aspects of structure, such as 0''" ground states for even-even 

^ I nuclei, and even vibrational spectra, arise with randomly distributed interactions. 

Ch ■ It is the purpose of this Letter to confront calculations with random and non-random 

interactions with robust experimental results to study which pervasive features of nuclei arise 

K> ! from the basis space and the generic type of interaction (e.g., 2-body) and which depend 

H [ on the details of the interaction. We will use the insights developed in this analysis to 

better understand structural evolution and its relation to shell structure, and to project the 

behavior of exotic nuclei. 

First, it is useful to succinctly summarize existing studies. In Shell model |I|-^ and 
IBM 0-^ approaches for many-body systems, the majority (typically ~70%) of calculated 
ground states have J^ = 0"*" even though 0^ states comprise a much smaller fraction of 
the basis space than other angular momenta. While it might be suspected that such a 
"pairing" property could emerge as a consequence of time reversal invariance, it has been 
shown [^ that this is not the case. Comparison of ground state wave functions in nuclei 
differing by two nucleons shows further evidence for a pairing relationship, reflecting the 
generalized seniority scheme |p. The analysis also extends to excited states. In the Shell 
Model, modest probabilities were found [|l[ for vibration-like spectra, defined by energy level 
ratios -R4/2 = E{4:f)/E{2f) near 2.0. However it was not possible to produce significant 
numbers of spectra with rotational character, R4/2 ~ 10/3. 

In the IBM analyses [i-IO] a similar preponderance of O''" ground states was found. Here, 



however, in contrast to the Shell Model, the interactions in an s-d boson space lead rapidly, 
with increasing boson number Nb, to large probabilities for both vibrational [-R4/2 ~ 2.0] and 
rotational [-R4/2 ~ 3.33] spectra. Indeed, the IBM treatments produced no concentrations 
of R4/2 values near any other values, as noted in Refs. @-§| . These results are reminiscent 



of empirical correlations @,|T0[ of 4j' and 2J' energies which give evidence for anharmonic 
vibrator and rotor behavior. 

These studies |I|-§[ are both surprising and important, and potentially impact our under- 
standing of the origins of some of the most characteristic features of nuclear spectra. Such 
features seem to emerge primarily from the nature of the basis space and the rank (2-body) 
of the interaction. 

Despite all this recent work, there has been no explicit comparison with experiment. It is 
therefore one of the purposes of the present Letter to do so, with data on level energies span- 
ning the entire nuclear chart ||Tl|]. This is the first time that detailed data on the low-lying 



nuclear spectra have been compared with predictions obtained from random ensembles. In 
the process, some previously unrecognized features of the data themselves will be discerned, 
which are not in fact produced with random interactions, giving us the kind of information 
needed to define favored interactions. 

To proceed, we use two IBM Hamiltonian ensembles, the full, most general, Hamiltonian 
(denoted Hjbm) with two 1-body and seven 2-body interactions, and the more "focused" 



Hamiltonian iJ^K = ^nd — kQ ■ Q of the extended consistent Q formalism (ECQF) |T2|,|T3 
where Q = s^rf + Ss + x[d^d]^'^^ ■ We use either uniform, random x ^ [0, — V^/2] or fixed 
X, as specified in each set of calculations below. The ensembles are defined by choosing 
random Gaussian interaction strengths in the Hamiltonians Hjbm and H^^ as in ref. [Q. 
The full Hamiltonian gives results essentially identical to those in Ref. [4], with peaks at 
-R4/2 ~ 2.0 and 3.33, and growing probability for the latter as A''^ increases. The focused 
ECQF Hamiltonian (with random x in the Q operator) gives similar but "cleaner" R4/2 
distributions. Both results are shown in the inset to Fig. 1. It is interesting to examine the 



correlations of -^(4^") vs. E{2i). Empirically, this correlation has been shown P,p!0[l to be 
extremely robust, exhibiting a bi-modal character with slopes of 2.0 and 3.33 in collective 
nuclei. Comparisons of these data with calculations using the full Hibm for boson number 
Ns=16 are shown in Fig. 1. The overall energy scale is fixed to the data by choosing the 
width of the Gaussian distributions. Similar results are obtained for other Nb- 

While the overall agreement is excellent, the trends of Fig. 1 only show part of the 
story. They do not easily reveal the density distribution along the trajectory, that is, the 
probability distribution of nuclei as a function of -R4/2- And it is from this perspective that 
specific empirical features appear that are not described by calculations unless the range of 
interactions is constrained. To see this, consider Fig. 2(left) which shows empirical i?4/2 
values for the entire nuclear chart from Z=8-100. In light nuclei, there is little evidence for 
collective R4/2 values. In the lowest medium mass region, 28 < Z < 50, however, a peak at 
R4/2 ~ 2.3 begins to emerge. In heavier nuclei this peak remains but is now accompanied 
by nuclei with larger i?4/2 values and, in particular, a sharp enhancement near R4/2 ~ 3.33. 
Fig 3a shows the composite distribution of -R4/2 values for all Nb- The peak at R4/2 ~ 2.3 
persists. 

Is this preference a fundamental feature of nuclear spectra or a result of some bias in 
the data? To study this we consider in Fig. 2 (middle) a different cut through the available 



i?4/2 data, namely, distributions for different ranges of boson numbers (i.e., valence nucleon 
numbers: N^ai = SiV^). Pre-collective nuclei {Nb < 3) show few i?4/2 values above 2.0. 
Nuclei with Nb > 13 are nearly all rotational, and those with Nb in the range of 4-12 are 
transitional. Surprisingly, no range of valence nucleon number shows an enhancement at 
-R4/2 ~ 2.0, contrary to the results of the random ensembles. (In this context we recall 
that the slope of 2.0 in Fig. 1 does not imply R4/2 = 2.0 due to the finite intercept of 
this segment.) However, the second panel of Fig. 2 (middle) shows that the peak centered 
on R4/2 = 2.3, mentioned above, is specifically associated with nuclei with 8 — 18 valence 
nucleons. To our knowledge, this striking feature has not been pointed out explicitly before 
but, as we shall see, leads to an empirical constraint on interactions that is not evident from 
Fig. 1 alone. 

Part of the explanation of the shell-by-shell results in Fig. 2(left) therefore merely reflects 
changing shell sizes and hence the distribution of possible Nb values. However, this is not 
sufficient. It does not, for example, explain the overwhelming preference for rotational i?4/2 
values in regions like 82 < Z < 102 where nuclei with small numbers of valence nucleons 
also exist, or the abundance of -R4/2 values near 2.3 in the 50 < Z < 82 shell. The data also 
reflect which nuclei are known in each region. For example, most known nuclei with Z > 82 
have large Nb- 

If we assume the universality of the empirical distributions in Fig. 2(middle), we can 
construct expected -R4/2 distributions for any region of nuclei simply by tallying the number 
of occurrences of a given A^^ for a range of N and Z values and multiplying that number 
by the relevant empirical -R4/2 distribution. In Fig 3b, we show the -R4/2 distribution for 
all nuclei with Z<82 assuming that, for each proton shell, all i?4/2 values are known for 
neutrons filling both the same shell and the next shell (regardless of where the drip lines 
are). So, for example, for 8 < Z < 20, we compile the -R4/2 distribution for 8 < A^ < 28. 
For 28 < Z < 50, we use 28 < A^ < 82. The peak at -R4/2 ~ 2.3 persists, reflecting simply 
the fact that boson numbers from 4-9 appear frequently in virtually all shells. Fig. 3c shows 
similar results for 82 < Z < 126 and 82 < A^ < 184. The peak at -R4/2 ~ 2.3 remains, but 
with a proliferation of larger boson numbers, a large abundance of rotational R4/2 values 
near 3.33 also appears. 

It is an important component of this Letter to see if calculations with random interactions 
can reproduce these -R4/2 distributions. Such calculations (for both Hjbm and H^,^) are shown 
in Fig. 2 (right) for comparison with the data in Fig. 2 (middle). Some results agree with 
the empirical trends such as the lack of structure for low boson numbers with Hjbm and 
the abundance of R4/2 ~3.33 values for large Nb- However, these random calculations tend 
to give a peak at -R4/2 ~2 and nowhere do they show a peak at 2.3. 

We will see momentarily that this points towards the need for speciflc non-random in- 
teractions. First, though, it is important to note that, although the interactions used in 
iJeK are random, they actually contain an implicit bias in R4/2 values. For Nb = 7 and e/n 
values from 00 down to ~ 50, IBM spectra are vibrational. For e/zt < 20 rotational spectra 
result. Thus, for a random interaction, which effectively samples an inflnitely large range 
of e/n values, the measures of the vibrational and rotational regions far exceed that of the 
transitional region, which accounts for the high frequencies of calculated R4/2 values near 2.0 
and 3.33. (Similar arguments apply, but for shifted ranges of e/n, for the other Nb values.) 
The Gaussian random numbers also enhance the 0(6) values of R4/2 (~ 2.5) for small Nb, 



as seen in Fig. 2 (right column, dashes). 

We carried out a large sampling of calculations to determine the specific range of inter- 
actions that yield a given i?4/2 value and, in particular, i?4/2 ~ 2.3 for Nb =4 — 9. Figure 4 
shows the statistical ranges of e/n and x values that give R4/2 in the ranges of 2.2-2.4 and 
2.8 - 3.0 for Nb=7. Given the discussion just above, it is not surprising that -R4/2 ~ 2.3 
values require intermediate ranges of e/n near 20-30. 

Thus, the empirical preference for i?4/2 ~ 2.3 values refiects two equally necessary in- 
gredients: the relevant interactions, and the locus of accessible nuclei. Had the valley of 
stability lain elsewhere, or had the e/n and x values been different, the i?4/2 distribution 
would have been different. Comparison of measured -R4/2 distributions in new, unmapped 
regions with expectations (as in Figs. 3(b)-(c)) can thus be used to deduce information 
either on shell structure (e.g., magic numbers) or on the interactions applicable to these new 
nuclei. 

Another example of characteristic nuclear data that constrains the interaction comes 
from rotational nuclei, namely, the energies of gamma vibrational modes, i?(2^), in deformed 
nuclei. In Fig. 5a we show the results for Hjbm with A^^ = 16 (to approximate the mean 
boson number for the deformed rare earth nuclei). There is no evidence for a peak in 
R2/2=E{22) / E{2f) at any value larger than 2. In heavy nuclei, however, such excitations 
are well known to occur at 10-20 times the 2^ energy. Fig. 5b shows a collection of empirical 
results for R2/2 for the rare earth region showing a broad bump centered on R2/2 ~ 15. If we 
now truncate the interaction space to that relevant for deformed nuclei, that is, e/n < 26 we 
see, as shown in Fig. 5c, that, for each x value, a distinct peak in R2/2 probability appears. 
The peak for x = —0.4 (solid), which is the traditional value for deformed nuclei, gives 
an excellent reproduction of the data. Once again, therefore, specific data on actual nuclei 
help identify and constrain the nature and relative importance of the 1- and 2-body residual 
interactions. 

To summarize, we presented a confrontation of the data for Z=8-100 with the predictions 
of nuclear structure calculations with random interactions, in the framework of the IBM. 
(a) General features of nuclear spectra, such as 0"*" ground states and the appearance of 
vibrational and rotational nuclei 0, and the global -^(4^^) vs. E{2f) trajectory, are typical 
of calculations with random interactions; (b) Heretofore unrecognized features of the data 
such as the absence of an abundance peak at the vibrator value i?4/2 = 2.0 and a frequency 
enhancement at -R4/2 ~ 2.3, in nuclei with intermediate numbers of valence nucleons, were 
discovered, (c) Such robust features of nuclear data, or others such as the characteristic 
energies of 7- vibrational modes in deformed nuclei, point to specific ranges of interactions. 
In particular, we showed that data such as these depend both on the relevant valence nucleon 
numbers (i.e., the nuclei that are known) and on the interactions. If the residual interactions 
change in new nuclei (e.g., near the drip lines), then so will the distributions of R4/2 values 
and the collective features of exotic nuclei. If the locus of newly accessible nuclei favors 
certain shell regions, the distribution of -R4/2 values and the nature of collective modes will 
refiect this as well. Conversely, valence nucleon number, and hence shell structure and the 
locations of magic nuclei, can be estimated from empirical knowledge of -R4/2 values, even 
when the magic nuclei themselves are not accessible. The usefulness of this result in studies 
of exotic nuclei, where issues of shell structure or the nature of residual interactions are 
paramount, is obvious. 
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FIG. 1. Plot of E{4f) against E{2f) for collective nuclei with 38 < Z < 82 (o), similar to 
[9]. In this and the following figures, we use all available data from [11]. The data lie along lines 
with slopes 2.0 and 3.33 (dots). The finite intercept of the former implies -£^(4^") = E{2~[) + e. 
Hence -R4/2 = 2 + e/E{2J) varies along this segment. The statistical distribution of E{4'{) and 
E{2i) energies for the full Hjbm is a strongly correlated function, whose maximum is indicated 
by the solid line, and the approximate FWHM by the dashes, which agrees well with the data. 
The variance of the Gaussian random numbers (which scales both axes) is adjusted to match the 
transition point between the two empirical slopes. (Inset) Statistical distribution of R4/2 for the 
full Hamiltonian (dots) and -ffe^ (solid) showing peaks at 2.0 and 3.33. (In the main figure, the 
points corresponding to the abundance peak at R4/2 ~ 2 are at high £^(2^), off'scale at the upper 
right.) All simulations are for Nb = 16 and 50000 realizations of the random Hamiltonians. 




FIG. 2. (Left column) Experimental distribution of R/^n as a function of proton shells. (Middle 
column) Experimental distribution of i?4/2 as a function of Nb- (The number of valence nucleons 
is 2Nb) (Right column) Statistical predictions of i?4/2 for selected boson numbers using Hjbm 
(solid) and H^^ (dashes) with random x- The vertical scale for this column is arbitrary. 
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FIG. 3. Empirical distributions of i?4/2 for indicated proton regions, (a) Sum of existing data 

from Fig. 2 (left), (b) and (c) Semi-empirical estimates of i?4/2 distributions for large regions of 

nuclei where it is assumed that the i?4/2 values would be available for full proton and neutron 

shells. Specifically, for each full proton shell, we take the set of nuclei obtained when the neutrons 

fill the same and next higher shell. These panels are constructed by assuming that the empirical 

distributions in Fig. 2 (middle) apply to all nuclei and, therefore, by folding these distributions 

with the frequency of Nb for these regions. 
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FIG. 4. Statistical regions of x ^nd ^/'^ where -R4/2 falls in the range 2.2-2.4 and 2.8-3.0 for 
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FIG. 5. (a) Distribution of R2/2 for Hjbm (dots) and H^,^ with random x (solid), (b) Experi- 
mental distribution of -R2/2- (c) Distributions for H^f^ for e/n < 26 and x = — v7/2 (dashes), -0.8 
(dots) and -0.4 (solid), showing the selectivity in i?2/2- 



